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Abstract 

We study a generalization of the standard trapping problem of random walk theory 
in which particles move suhdiffusively on a one-dimensional lattice. We consider 
the cases in which the lattice is filled with a one-sided and a two-sided random 
distribution of static absorbing traps with concentration c. The survival probability 
that the random walker is not trapped by time t is obtained exactly in both 
versions of the problem through a fractional diffusion approach. Comparison with 
simulation results is made. 
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1 Introduction 



The trapping of Brownian particles by static traps randomly distributed over 
either a Euclidean or a disordered substrate is a fundamental problem of non- 
equilibrium statistical mechanics and chemistry with a very wide range of 
applications [1,2,3,4,5]. This is also one of the oldest problems in random walk 
theory which was essentially formulated by Smoluchowski at the beginning of 
past century in his theory of coagulation of colloidal particles [1,2,3,5]. This 
model has proven useful in research areas such as the trapping of mobile de- 
fects in crystals with point sinks [6,7,8], the kinetics of luminescent organic 
materials [7], the kinetics of photosynthetic light energy to oxygen conver- 
sion [9], anchoring of polymers by chemically active sites [10], atomic diffusion 
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in glasslike materials [11], and others [12]. Some generalizations of the ba- 
sic model have also been considered recently. The case of a gated (ungated) 
random walker trapped by a distribution of ungated (gated) fixed traps has 
been proposed in connection with the kinetics of reactions between complex 
molecules whose active groups are screened as the molecules diffuse [13]. The 
trapping problem with many random walkers has also been studied in one 
dimension [14] and the corresponding trapping statistics have been found on 
Euclidean and fractal lattices [15,16,17]. The multiparticle predator-prey prob- 
lems in which a single fixed trap, the "lamb" , is captured by one of a set of N 
random walkers or "lions" initially placed at a given distance from the prey 
was also discussed by Krapivsky and Redner [18,19]. Another quite interest- 
ing and difficult variation of the standard trapping problem is that where the 
traps diffuse too [5,20]. 

Trapping reactions between molecules embedded in biological samples and 
disordered materials are usually handicapped by the porous and statistical 
fractal structure of these media [4]. In some cases this gives rise to subdiffusion 
of the particles, i.e., the mean square displacement (r^(t)) of the particles from 
the original starting site is no longer linear on time, but verifies a generalized 
Pick's second law: 

where 7 (with < 7 < 1) is the (anomalous) diffusion exponent and is 
the diffusion coefficient. Of course, there are many other instances in which 
subdiffusion processes appear [21,22,23,24,25]. A useful approach for under- 
standing subdiffusion processes is by means of the continuous time random 
walk model (CTRW) in which the random walker performs jumps with a wait- 
ing time distribution with a broad long-time tail: ^(t) ~ ^-(i+t) for large t 
[1,2,3,21,26,27]. The long-time tail of this waiting time distribution incorpo- 
rates, in a statistical sense, the effect of the bottlenecks and the dead-ends in 
the diffusion of the random walker embedded in the disordered structure, and 
the model is compatible with Eq. (1). For subdiffusive random walkers the 
continuum description given by the ordinary diffusion equation is replaced by 
the fractional diffusion equation [27] 

^W{x,t)^K,oDl-'^^W{x,t), (2) 

where oDI~'^ is the Riemann-Liouville fractional derivative of order 1 — 7 
[27,28,29], 

Dl-m{x, t) = -i-l- /* ^^ ^(^^j) 

* ^ ' ^ r{j)dtJo {t-ry-^' ^ ^ 

and W{x, t) is the probability density that the particle that started at at 
time is at a; at time t. 

Our objective in this paper is to find analytical expressions for the survival 
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probability defined as the probability that no trap site has been reached 
by the suhdiffusive random walker by time t. Two variations of the classi- 
cal trapping problem (sometimes called Rosenstock's trapping problem) are 
considered: (i) the "one-sided" trapping problem [14] in which only one half- 
line of a one-dimensional lattice is filled with a random distribution of static 
traps (this process could mimic the excitation or production of defects on one 
side of a fiber by irradiation, the other side being shielded); and (ii) the "two- 
sided" trapping problem corresponding to the trapping of a single subdiffusive 
random walker placed initially at x = between two half-lines (x < and 
X > 0) randomly filled with static traps. We find exact analytical solutions for 
both problems and compare them with simulations. As precedents of these re- 
sults we whould cite the work of Blumen, Klafter and Zumofen (see [12,30,31] 
and references therein) about the trapping of particles on Euclidean and frac- 
tal substrates for random walkers with waiting time distributions with broad 
long-time tails: '^{t) ~ t"^"^ , < 7 < 1. 

The paper is organized as follows. In Sec. 2 the exact solution of the one- 
sided trapping problem is obtained and compared with simulation results. 
In Sec. 3 we build the extended Rosenstock's approximation for the survival 
probability of the particle in the one-sided trapping problem by calculating 
an exact expression for the moments of the one-sided span (the number of 
distinct sites visited by the random walker in a given direction). In Sec. 4 we 
obtain the exact survival probability of the subdiffusive random walker for the 
two-sided case in integral form. This integral is evaluated in the asymptotic 
long-time and short-time limits in order to get simple closed expressions. The 
paper ends with some conclusions and remarks in Sec. 5. 



2 The one-sided trapping problem 



In the one-sided trapping model quenched traps are randomly distributed on, 
say, the right-hand side of a one-dimensional lattice {x > 0) with concentration 
c. The random walker is placed initially upon site x — Q and performs jumps 
to its nearest neighbour sites with a waiting time distribution ■0(t) until it 
reaches a trap site where it is absorbed. The survival probability is given by 

[1]: 

m-i:e-^'P{t\r), (4) 

r=l 

where A = — ln(l — c) and P{t\r) is the probability that the span of the 
random walker in the positive direction (the largest distance reached by the 
random walker for a; > 0) is equal to r after t time steps. Now let T{t\r) be 
the probability that the site x — r has not been visited by the random walker 
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by time t (the so-called fixed-trap survival probability) . Because [32] 

dT{t\r) 



P{t\r) 



dr 



(5) 



Eq. (4) takes the following form in the continuous limit: 

m 

or, integrating by parts, 



dr 



(6) 



<5(t) = A / e-^''r(t\r)dr , 
Jo 



(7) 



because r(i|0) = 0, r(f|oo) = 1 and e — > as r ^ oo. For subdiffusive 
particles, the function T{t\r) can be written in terms of Fox's H function [33]: 



r(t|r) = l-H 



r 


(1,7/2)' 




(0,1) 



(8) 



The time Laplace transform of r(t|r) is given by 



1 



rfslr) = - 1 — exp(— r\/s^) 

S L 



(9) 



Then the Laplace transform of ^{t) is readily calculated from Eqs. (7) and 
(9): 

^(«) = — r4^i—i (10) 



s + X^K-yS^-^/^ 
so that [34] 

m = E,/,{-0 , (11) 

where ^ = X^Jlc^ and Ea{z) is the Mittag-Leffler function with parameter 
a [34,28]. For 7 = 1, the Mitagg-Leffler function becomes [27,28]: 

$(i) = e«Wc(e) , (12) 

with ^ = \\fDi, {Ki = D), and we recover the result for the survival probabil- 
ity of a normal diffusive random walker in the one-sided trapping problem [14]. 
For very long times, the asymptotic expansion of the Mittag-Leffler function 
[28] allows us to write 



^^'^ ,tir(i-^7/2)^ ^ 



(13) 
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Fig. 1. Survival probability versus ^ = Xy^Kyfy. The circles are simulation 
results for subdiffusive walkers with 7 = 1/2. The squares are simulation results 
for normal diffusive walkers with ip(t) = e~*. The lines are the corresponding exact 
results given by Eqs. (11) (solid line) and (12) (dashed line), respectively. 



Thus, an asymptotic time regime is reached for t ^ l/{Ky\'^Y^'^ where the 
survival probability exhibits a power law decay 



r(l-7/2)A/^' 



(14) 



This is an a/^feftraic fluctuation slowdown corresponding to the Donsker-Varadhan 
limit [1,2,35]. We compare in Fig. 1 the exact survival probability for a subd- 
iffusive random walker with a waiting time distribution with the form of the 
Pareto law 
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(15) 



and 7 = 1/2. Taking the distance between nearest neighbour sites (jump 
length) as 1, the subdiffusion constant is K1/2 = 1/ [2r(l/2)] = I/a/Itt. The 
concentration of traps used was c = 0.01. 



3 The Rosenstock's approximation for the one-sided trapping prob- 
lem 



Let S{t) be the number of distinct sites on the positive half-hne visited up to 
time t by a random walker who started at a; = at time t = 0. Then, the 
survival probability of the random walker is given by = ((1 — c)"'''^*^^ = 
^g-A5(t)^^ the average being performed over all realizations of the random 
walker's exploration of the lattice from time until time t. By means of the 
cumulant expansion technique [1,2] the following equivalent form is derived: 



^{t) = exp 



E(-ir 



.n=l 



(16) 



where Kr 



n 



1,2,... denote the cumulants of S{t): Ki = {S{t)), K2 



(S"^ (t)) — {S (t)) = cT^(t), .... If we keep the first n+1 terms of the sum in equa- 
tion (16) we arrive at the nth-order Rosenstock approximation [1,2,30]. The 
error made by using this approximation is 0{Kn+2^"'~^'^)- Thus, the condition 
K'n+2 <^ 1/ ^"'^'^ must be fulfilled for the nth-order Rosenstock approximation 
to be reasonable. 



For the Rosenstock approximation to be useful for the subdiffusive case, it 
is necessary to know the cumulants or, equivalently, the moments of S{t) for 
subdiffusive particles. The next step is thus to derive exact expressions for the 
moments (S'™(t)), m = 1, 2, ... of the one-sided span of a subdiffusive random 
walker, i.e., the territory explored in a given direction. The probability density 
of the one-sided span, P{t\r), was given in Eq. (5) in terms of the survival 
probability. Hence 

{S^{t)) = r ^^^^r'^dr ^ - r ^(1 - r(t\r))r"'dr (17) 
Jo dr Jo dr 

or, integrating by parts, 

(S"'(t)) = m / [1 - r(t\r)] r'^-^dr , (18) 
Jo 

where we have taken into account that r(t|0) = and limr^^r"^[l — T{t\r)] — 
as a consequence of the stretched exponential behaviour of 1 — r(t|r) as 
r — > 00 for subdiffusive particles [33]. Prom Eq. (18) and the exphcit expression 
of the survival probability in Eq. (8) we obtain 

{S^{t))=ms^{K,t''r/' , (19) 

with 

dz . (20) 



^m-l 0-10 







z 






(0,1) _ 
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In order to calculate the coefficients Sm, m — 1,2,... we observe that the 
Laplace transform of the Fox H function appearing in the integrand in Eq. 
(20), H{u), is a generating function of these coefficients: 



u 



m—l 



(21) 



Taking into account some properties of the Fox functions [36,37,27] we finally 
identify H{u) with a two-parameter Mittag-Leffler function as follows: 



H{u) = £;^/2,l+^/2(-1i) , 

where, in the last identity, we have used the relation [27] 



(22) 



■al3[ 



-u) = H 



U 



(0,1) 

(0,1), (1-/3, a) 



(23) 



between the two-parameter Mittag-Leffler functions and the Fox functions 
[27]. The function Eap{u) admits the following series expansion [27,28]: 



V{ak + /5) 
From Eqs. (21), (22) and (24) we find 

_ (m-1)! 
r(l + m7/2) 

so that from Eq. (19) we finally obtain 

m! 



a > 0, /3 > . 



(24) 



m = 1,2,... 



(25) 



r(l + m7/2) 



m/2 



(26) 



Obviously, for 7 = 1 we recover the exact result for the moments of the one- 
sided span of a diffusive random walker [14]: 



(5-(t)) = r[(^+l)/2] 



TT 



(27) 



From Eq. (26) one can calculate the cumulants Kn,n — 1,2,. . . of the one-sided 
span distribution, and their direct substitution into Eq. (16) yields the general 
expression for the Rosenstock approximation for subdiffusive particles: 



m = exp 5: ^—^anC 



.n=l 



(28) 
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The first three coefficients are: 



r(i + i) 

1 2 

'''~~[r(i + i)]2 + r(i + 7) ^^^^ 

_ 2 6 6 

" [r(i + i)P " r(i + i)r(i + 7) + r(i + ^) ' 

It is interesting to note that from Eqs. (11) and (28) we get the following 
identity 

ln£;,/,(0 = f:^r, (30) 

n=l 

which provides a series expansion of the logarithm of the Mittag-Leffler func- 
tion. 



4 The two-sided trapping problem 



In this section we calculate the exact survival probability of a single subdif- 
fusive random walker placed initially between two half-lines populated with a 
random distribution of traps with concentration c. We will use parallel argu- 
ments to the standard ones for diffusive random walkers [2,38]. We must also 
mention that an alternative approach was proposed by Anlauf [39] to find the 
long-time behavior for the one-dimensional trapping problem with a diffusive 
particle. 

As a starting point we evaluate the probability W{x, t\xo, t = 0) that a random 
walker starting from xq att — inside a box < a:; < L with absorbing bound- 
aries is at X at time t. The function W{x,t\xo,t = 0) satisfies the fractional 
partial differential equation (the subdiffusion equation): 

0'^ 

-W{x, t\xo, Q) = K^o Dl-"'—W{x, t\xo, 0) , (31) 

with the boundary conditions W{Q,t\xo,0) = W{L,t\xo,0) = 0, and the ini- 
tial condition W{x,0\xo,0) — 5{x — Xo) with < 7 < 1. Equation (31) is 
straightforwardly solved by separation of variables [27] : 

W{x, i|xo, 0) = ^ E sin (-^ j sin (^-^ j \-K^—-tA , (32) 

where we have taken into account that the solution of the fractional differ- 
ential equation dT/dt = —K^}? qDI^^T is given by a Mittag-Leffler func- 
tion T{t) = {—K^X'^P). The survival probability of a subdiffusive random 
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walker starting from xq G (0, L) is given by Jq W(x,t\xo,0)dx. Let 
be the survival probability of the random walker averaged over all config- 
urations of traps which contain the origin inside a hole of length L. Then 
= J^W{x,t\xo,0)dxdxo, and Eq. (32) yields 

(33) 

The survival probability of the subdiffusive random walker in the two-sided 
trapping problem is finally given by an average of over the distribution 

of hole lengths in the random trap configurations. This distribution is rj{L) — 
\^Le~^^ [2], and consequently we have 

/■oo 

= AM Le-^^^L{t)dL 

with f3n = K^{2n + l)^7r^t^, n = 1, 2, .... This is the exact integral represen- 
tation of the survival probability of the subdiffusive particle in the two-sided 
trapping problem. 



4-1 Long-time behavior of ^{t) 



For 7=1 the Mittag-Leffler function reduces to an exponential and the inte- 
grals in Eq. (34) arc estimated asymptotically for large times by means of the 
Laplace method [2]. However, this approach is not possible for 7 7^ 1. In order 
to perform a long-time asymptotic evaluation of the integrals in Eq. (34), we 
use the following series expansion of the Mittag-Leffler function: 

which is valid for I3n/L'^ — > oo in the range < 7 < 2 [28]. Inserting this series 
expansion into Eq. (34), integrating term by term, and after some rearrange- 
ments we find 

m) - - f (4-2-^"-)C(2m + 2)r(2m + 2) . 




m=l 



where we have taken into account that E^=o l/(2n+l)^"'+2 = (l-2-2"^-2)C(2m+ 
2), C(s) being Riemann's zeta function. The dominant term is 
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which is vahd for ,^ ^ 1 . Therefore we find an algebraic fluctuation slowdown 
of the survival probability corresponding to the Donskcr-Varadhan [35] limit 
with an exponent two times greater than that of the one-sided long-time be- 
havior in Eq. (14). The two-sided result ~ t~'^ is well known and has 
been interpreted as an "avoided crossing effect" : trapping cannot be more effi- 
cient that the probability ilj{t)dt ~ t~'^ of remaining at the initial site (see 
[12,31,30] and references therein). However, note that this "avoided crossing 
effect" is absent for the one-sided trapping problem since the probability of 
trapping <l>(t) ~ f^/^ decays even more slowly than the probability of remain- 
ing at the initial site. 



4-2 Short-time behavior of ^{t) 



A simple analytical expression for the short-time behavior of <l>(t) can also 
be derived. To do so, we calculate the time derivative of $l(^) in Eq. (33) 
obtaining 

|$L(t) = -^,K, oA'-V(7) , (38) 

where we have taken into account that Mittag-Leffler functions are the solution 
of the fractional relaxation equation [27] and 



oo 

J(7) = E^7(-4a'(n+l/2)- 



n=0 



(39) 



with a^ — K^n'^V / L?' . The major contribution to the sum defining 7(7) in Eq. 
(39) when a is small comes from large n. Hence, to find the lowest-order term 
in the small-a behavior of ^(7), we can approximate the sum by an integral: 



J(7) - / EJ-Aa^x^]dx = —EJu = 0) , 
^0 2a 



(40) 



where E^{u) = e '^yE^{y^)dy is the Laplace transform of E^{y^). To cal- 
culate E^{u) we will exploit the relation between the Mittag-Leffler and the 
Fox H functions [27,37]: 



(0, 1/2) 

(0,1/2), (0,7/2) 

Prom Eq. (41) and the properties of Fox functions [36,37] we find 



(41) 



21 



-^7 ("5) = 2^22 



(1/2, 1/2), (1-7/2,7/2) 
(0, 1), (1/2, 1/2) 



(42) 
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Inserting the result in Eq. (42) into Eq. (40), we finally arrive at the following 
explicit expression for J(7): 



where we have calculated the value of the Laplace transform for s = by 
resorting to the general series expansion of Fox functions [27]. From Eqs. (38) 
and (43) we can also write 

dt ' r(7/2) L ' ^ ' 

where the fractional derivative has been performed according to the rule: 
oD^ = r(l + /i) t^~'^/r(l +/U — z/) which is valid for arbitrary real parameters 
fj, and z/ [27]. By integrating Eq. (44) with the initial condition ^L{t — 0) — 1 
we find 
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and performing the average in Eq. (34) over trap configurations we finally 
obtain the short-time behavior of the survival probability of a subdiffusive 
random walker in the two-sided trapping problem 

or, for i <^ I, 

*W~-p(-r(T^)- ^"'> 

This expression is just the zeroth-order Rosenstock approximation $(i) = 
exp[— A(S'^(t))] [see Eq. (16)] because the average two-sided span {S^{t)) 
is simply twice the average one-sided span {S(t)) and, by Eq. (26), {S(t)) = 
{K^py^'^ /T{l+'y /2). In Fig. (2) we compare simulation results for subdiffusive 
random walkers in the case 7 = 1/2 with the numerical integration of the exact 
integral representation of ^{t) given by Eq. (34), the long-time behavior given 
in Eq. (36), and the Rosentock short-time behavior in Eq. (47). 



5 CONCLUSIONS AND REMARKS 



In this paper we considered the one-dimensional trapping problem for subdif- 
fusive particles. In this problem a single subdiffusive random walker starting 
at a; = on a lattice randomly filled with absorbing static traps with a den- 
sity c performs a random exploration until he encounters one of these traps. 
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Fig. 2. Survival probability <I>(t) versus ^ = Xy^K^tf for subdiffusive random walkers 
with 7=1/2 for the two-sided trapping case. The symbols are simulation results for 
c = 0.01 and 5 x 10^ trials. The solid line corresponds to the numerical integration 
of the exact expression (34). The dashed and dash-dotted line are the long-time 
behavior predictions obtained by keeping a single term and three terms in Eq. (36), 
respectively. The dotted line is the short-time behavior given by Eq. (47). 

The quantity of main interest is the survival probability (f){t) of the particle. 
We have considered two versions of the problem: the one-sided case in which 
only one half-line is filled with traps, and the two-sided case corresponding 
to a random filling of both sides of the line x > and x < 0. In the context 
of ordinary diffusive random walkers, this problem has a long tradition, and 
many applications to physics and chemistry have been discussed [1,2,3,4]. 

Great interest has also arisen recently around subdiffusive anomalous processes 
(for example, as a way of mimicking transport in disordered media), and a 
continuous fractional diffusion description has been put forward (see [27] and 
references therein). It thus seems convenient to extend the fruitful trapping 
model to the case of subdiffusive particles. We achieved this objective for the 
trapping of a single subdiffusive random walker in one dimension and, by 
means of the fractional diffusion formalism, exact analytical expressions were 
found in terms of the special functions characteristic of fractional calculus. 
There is a possibility for these results to be checked in diffusion experiments 
performed in constrained geometries and/or disordered materials [25]. From 
another point of view, the trapping problem could also be interpreted as a 
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pseudo-first-order reaction of the form A + T ^ T, wfierc A is tfie random 
walker and T is tlie trap. Tliere is also an increasing interest in reaction- 
subdiffusion processes [40,41], and our exact results apply to a special class of 
these processes. 



The long-time behavior of the survival probability has been a subject 
of particular interest in the ordinary diffusive trapping problem [1,2,35,39]. 
In this regime an anomalous fluctuation slowdown, known as the Donsker- 
Varadhan limit, has been reported. In the one-sided and two-sided subdif- 
fusive models we found the algebraic decays <I>(t) ~ t^^/^ and <I>(t) ~ t"''', 
respectively. For normal diffusive particles, the algebraic decay ~ 
was found in [14] for the one-sided trapping case and interpreted as due to the 
long explorations carried by the random walkers on the half-line free of traps 
[14]. A stretched exponential decay is the corresponding behaviour for the 
long-time two-sided trapping of diffusive particles [39]. The two-sided subdif- 
fusive result ~ is well known and has been interpreted as an "avoided 
crossing effect" (see [12,31,30] and references therein). This effect is absent for 
the one-sided trapping problem since the probability of trapping ~ t~'^/'^ 
for this case decays even more slowly than the probability of remaining at the 
initial site. 



The present work can continue along several directions. First, the generaliza- 
tion of the subdifTusive trapping problem to the case of > 1 independent 
subdiffusive random walkers could be analyzed as has already been done in the 
diffusive case [14]. Application of the same techniques to higher dimensional 
spaces is also an obvious approach, but we consider it unlikely to produce 
exact closed expressions for the relevant quantities. Trapping models of subd- 
iffusive particles in which the trapping process is stochastically "gated" is also 
an interesting field. The effect of "gates" in diffusion-limited reactions has re- 
cently been considered in connection with reactions of complex molecules in 
biological media [13]. The chemically active groups of these molecules may be 
screened by the inactive parts, giving rise to effective reactivities described by 
Poisson processes. The media where these reactions takes place are usually dis- 
ordered and rcaction-subdiffusion models should provide a better description 
of these processes. 
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